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A relativistic distorted wave impulse approximation model for electron-induced one proton 
knockout reactions and a Green's function approach to inclusive scattering are developed. 
Results for (e, e'p) and (e, e') reactions are presented in various kinematical conditions and 
compared (when possible) with nonrelativistic calculations. 
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>'. 1. INTRODUCTION 

O 

• A long series of high-precision experiments on several nuclei [1-7] have generated a well es- 

tabhshed tradition which singles out exclusive (e, e'p) knockout reactions as the primary tool to 
CN , explore the single-particle aspects of the nucleus. Theoretical calculations were carried out within 

the framework of a nonrelativistic distorted wave impulse approximation (DWIA), where final-state 
interactions (FSI) and Coulomb distortion of the electron wave functions are taken into account [8] . 
Phenomenological ingredients were used to compute bound and scattering states. This approach 
was able to describe to a high degree of accuracy the shape of the experimental momentum dis- 
, tribution for several nuclei in a wide range of different kinematics [9, 10]. However, a systematic 

rescaling of the normalization of the bound state, interpreted as the spectroscopic factor for the 
corresponding level, had to be applied in order to reproduce the magnitude of the experimental 
■ distribution. 

' Similar models based on a fully relativistic DWIA (RDWIA) framework were developed in more 

recent years. The Dirac equation is solved directly for the nucleon bound and scattering states [11- 
13] or, equivalently, a Schrodinger-like equation is solved and the spinor distortion by the Dirac 
scalar and vector potentials is incorporated in an effective current operator in the so-called effective 
Pauli reduction [14-16]. 

A successful description of new (e, e'p) data at higher momentum transfer from Jefferson Lab- 
oratory (JLab) [17, 18] has been achieved within the framework of RDWIA, but slightly differ- 
ent spectroscopic factors are deduced, because the relativistic optical potentials in general give a 
stronger absorption than the corresponding nonrelativistic ones [13, 19]. Moreover, the limits of 
validity of the nonrelativistic DWIA analysis versus RDWIA were not always properly explored, as 
' discussed in Ref. [16], resulting, e.g., in a certain degree of ambiguity for the spectroscopic factors 

extracted at low energy. 

In Sec. 2 a RDWiA approach to low- and high-energy (e, e'p) data is presented and a careful 
analysis of the limits of the nonrelativistic DWIA is carried out. The sensitivity to different off-shell 
prescriptions for the electromagnetic current operator will also be discussed [20, 21]. 

In Sec. 3 a Green's function approach for the inclusive electron scattering is developed [22]. 
In the inclusive (e, e') scattering only the scattered electron is detected whereas the final nuclear 
state is not determined. The one-body mechanism is assumed to give the main contribution to the 
reaction in the quasielastic region. However, when the experimental data of the separation between 
the longitudinal and transverse responses became available it was clear that a more complicated 
framework than the single-particle model coupled to one-nucleon knockout was necessary. A review 
till 1995 of the experimental data and their possible explanations is given in Ref. [9]. Thereafter, 
only a few experimental papers were published [23, 24]. New experiments with high experimental 
resolution are planned at JLab [25] in order to extract the response functions. 

From the theoretical side, a wide literature was produced in order to explain the main problems 
raised by the separation, i.e., the lack of strength in the longitudinal response and the excess of 
strength in the transverse one. The more recent papers are mainly concerned with the contribution 
to the inclusive cross section of meson exchange currents and isobar excitations [26-28], with the 
effect of correlations [29, 30], and the use of a relativistic framework in the calculations [28]. 



X 
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At present, however, the experimental data are not yet completely understood. A possible 
solution could be the combined effect of two-body currents and tensor correlations [29, 31, 32]. 

In our Green's function approach to (e, e') scattering the spectral representation of the single 
particle Green function, based on a biorthogonal expansion in terms of the eigenfunctions of the 
nonhermitian optical potential, allows one to perform explicit calculations and to treat final state 
interactions consistently in the inclusive and in the exclusive reactions. 



2. THE (e, e'p) REACTION 
A. Relativistic current 



The main ingredient of the calculation is the nuclear transition amplitude, i.e., 

= J drx,^~\r)j^ exp{iq ■ r}(p{r) . (1) 

In RDWIA it is calculated using relativistic wave functions for initial and final states. 

The bound state wave function, ip, is given by the Dirac-Hartree solution of a relativistic La- 
grangian containing scalar and vector potentials deduced in the context of a relativistic mean field 
theory that satisfactorily reproduces single-particle properties of several spherical and deformed 
nuclei [33]. 

The ejectile wave function, x^~\ is written in terms of its positive energy component ^ f-\- 
foUowing the direct Pauli reduction method 

X^-^=f -V ^] , (2) 

V M + E' + S-V ) 

where S = S{r) and V = V{r) are the scalar and vector potentials for the nucleon with energy 
E' [34] . The upper component x+ is related to a Schrodinger equivalent wave function by the 
Darwin factor D{r), i.e., 

,,^^),,, Oir)^^i±^^. (3) 

$ J is a two-component wave function which is solution of a Schrodinger equation containing equiv- 
alent central and spin-orbit potentials obtained from the scalar and vector potentials. 

The choice of the electromagnetic operator is a longstanding problem. Here we discuss the three 
current conserving expressions [15, 35] 

= ^i(Q')^ + 4^GM{Q')a^''q. , (4) 

where q'^ = {uJ,q) is the four-momentum transfer, =\ q p — w^, = {E + E',Pm n 
is the anomalous part of the magnetic moment, Fi and F2 are the Dirac and Pauli nucleon form 
factors, Gm = Fi + is the Sachs nucleon magnetic form factor, and a^"^ = (i/2) [y^,j'^]. These 
expressions are equivalent for on-shell particles thanks to Gordon identity. However, since nucleons 
in the nucleus are off-shell we expect that these formulas should give different results. Current 
conservation is restored by replacing the longitudinal current and the bound nucleon energy by [35] 

J'' = r = r^ j\ E= VI |2 +M2 = ^\p'-q |2 +M2 . (5) 
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B. Nonrelativistic current 



In nonrelativistic DWIA the transition amplitude of Eq. (1) is evaluated using eigenfunctions 
of a Schrodinger equation for both the bound and scattering states. In standard DWIA analyses 
phenomenological ingredients are usually adopted. In this work and in order to perform a consistent 
comparison with RDWIA calculations, we employ for the bound state the upper component of the 
Dirac wave function ip and for the final state the Schrodinger-like wave function ^j. 

The nuclear current operator is obtained from the Foldy-Wouthuysen reduction of the free- 
nucleon Dirac current through an expansion in a power series of 1/M. Performing the expansion 
through second order we get 



i(o) - ^1 ' - u , Jp) - ' 

7rm = , = la X q-\ F^P , i/n-, = — -^^ — ;r — - tua x P . (6) 



(1) - ' J(2) -i(T-Pxq) , 



C. The (e, e'p) cross section 



The coincidence cross section of the (e, e'p) reaction can be written as the contraction between 
the lepton tensor, completely determined by quantum electrodynamics, and the hadron tensor, 
whose components are given by suitable bilinear combinations of the nuclear transition amplitude 
in Eq. (1). In case of unpolarized reactions the cross section can be written in terms of four 
response functions, Rw, as 



a = K {vlRl + vtRt + vltRlt cos (■(?) + vttRtt cos (2'i?)} 



(7) 



where is a kinematic factor, and i? is the out-of-planc angle between the electron scattering plane 
and the {q,p') plane. The coefficients vw are obtained from the lepton tensor components and 
depend only upon the electron kinematics [9, 10]. The response functions are defined as 



Rlo,{J^ (jO)^) 



Rlt oc -2 Re 



(J-(j°)t) 



Rt oc {r ( + {jy (jy)^) , 
Rtt {r {r)^) - {jy (jy)^) , 



(8) 



where average over the initial and sum over the final states is performed fulfilling energy conserva- 
tion. In our frame of reference the z axis is along q, and the y axis is parallel to q x p'. 

If the electron beam is longitudinally polarized with helicity h, the coincidence cross section for 
a knocked out nucleon with spin directed along s can be written as 



1 

2"" 



1 + P ■ s + h(^A + P' ■ 



(9) 



where a is the unpolarized cross section of Eq. (7), P the induced polarization, A the electron 
analyzing power and P' the polarization transfer coefficient. We choose for the polarimeter the 
three perpendicular directions: L parallel to p', N along q x p', and T = N x L. 



D. Results and Discussion 

In Fig. 1 the unpolarized ^^0{e,e'p) reaction leading to the pi/2 ground state and the p^/2 

first excited state of ^^N is considered. In the left panel, data have been collected at NIKHEF in 
parallel kinematics at a constant proton energy of 90 MeV in the center-of-mass system [7] . They 
are presented in the form of the reduced cross section [9]. The results for the transition to the 
Pi/2 ground state have been multiplied by 40. The dot-dashed lines refer to the nonrelativistic 
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FIG. 1: Left panel: reduced cross section for the ^^0(e, e'pY^N transition to the pi/2 ground state and 
first excited state of ^'""N at Ep = 90 McV constant proton energy in the contcr-of-mass system in parallel 
kinematics [7]. Right panel: cross section for the same reaction but at = 0.8 (GeV/c)^ in constant {q,oj) 
kinematics [17]. Data for the pi/2 state have been multiplied by 40. Dashed, solid, and dotted lines represent 
the results of the RDWIA approach with ccl, cc2, cc3 off-shell prescriptions, respectively. Dot-dashed lines 
in the left panel are the nonrelativistic results. 



calculations of Sec. 2 B. The theoretical results have been rescaled in order to reproduce the data 
by applying the spectroscopic factors Zp.^^.^ = 0.64 and Zp.^^^ = 0.54, respectively. The solid lines 

show the results of the RDWIA analysis with the cc2 off-shell prescription; dashed and dotted lines 
indicate the results when using the ccl and cc3 recipes, respectively. The resulting spectroscopic 
factors, Zp^^^ = 0.708 and Zp^^^ = 0.602, have been obtained by a fit using the cc3 current 

21], which gives an overall better description of the {e. e'p) obscrvablcs. Only small differences are 
:'ound between the rclativistic and the nonrelativistic models. Thus, they are almost equivalent in 
comparison with the data, which are reasonably described by both calculations. 

In the right panel, the same reaction is considered at the JLab constant (g, lo) kinematics with 
= 0.8 (GeV/c)^ [17]. The data now refer to the differential unpolarized cross section. The 
Pi/2 results arc multiplied by a factor 40. The theoretical curves have the same meaning as in the 
left panel and are rescaled again by the same spectroscopic factors. Only RDWIA calculations are 
shown since at the proton energy of this experiment relativistic effects are large and a nonrelativistic 
analysis gives unreliable results. The agreement with the data is very good also in this case. This 
outcome is particularly welcome, since the spectroscopic factors correspond to a nuclear property 
that must be independent of Q^. 

In Fig. 2 the response functions measured at JLab in the same kinematics with = 0.8 
(GeV/c)^ [17] are displayed and compared with our RDWIA calculations. The agreement with the 
data is satisfactory and of about the same quality as in other relativistic analyses [13, 15], but for 
the Rlt response function, where only the cc3 calculation reproduces the pi/2 data at low missing 
momentum while the cc2 one better reproduces the data. In Fig. 3 the polarization transfer 
components P'^,P'^ are shown as functions of the missing momentum for the ^^0{e,e'p) 
reaction at = 0.8 (GeV/c)^ and constant (g,a;) for the transitions to the -"^^N Pz/2-, f'^d 
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FIG. 2: Response functions for the ^^0(e,e'p) reaction at — 0.8 (GeV/c)^ in constant (q,a;) kinematics 
[18] leading to the ^^N (left column) and P3/2 (right column) residual states. Line convention as in Fig. 



respectively [18]. For these observables and in this kinematics, the sensitivity to ofF-shell 

effects is at most < 15%. The overall agreement with the data is still good. 

As a first step to study the role of meson exchange currents (MEC), we have considered in Ref. 
[36] the contribution due to the seagull diagram. The seagull current is written in momentum space 
as in Refs. [37, 38] and with the cutoff A = 1250 MeV in the pion propagator. The inclusion of 
the seagull diagram enhances the RDWIA results, but, in contrast to (7,^) scattering, the effects 
are generally small and visible only at high missing momenta. Thus, the comparison with data, 
that were already well reproduced by the direct knoclout model, is practically unaffected. In 
particular, no significant effects were obtained on the polarization observables from MIT-Bates on 
i2C(e,e'p)^^B [39] and from JLab on ''^Oie.e' v f^^ [18] reactions. 



3. THE (e, e ) REACTION 

In the one photon exchange approximation the inclusive cross section for the quasielastic (e, e') 
scattering on a nucleus is given in terms of two response function as [9] 

ainc = K{2eLRL + RT) , (10) 
where .fC is a kinematical factor and 

£L = ^ (l + 2^tan2(^9,/2)j (11) 

measures the polarization of the virtual photon. In Eq. (11) "i^g is the scattering angle of the 
electron. The longitudinal and transverse response functions, Rl and Rt, contain all nuclear 
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FIG. 3: Polarization transfer components P'^,P'^ for the ^^0{e,e'p) reaction at = 0.8 (GeV/c)^ in 
constant {q,to) kinematics [18] leading to the ^^N pi/2, P3/2 and si/2 residual states. Line convention as in 
Fig. 1. 



structure information and are defined as in Eq. (8) . They are directly related to the hadron tensor 
components 

W''''iu;,q) =P\{^f\ J%{q) \ ^0)? 6{E^ + u - Ef) = \J';^ {q)G{Ef)J''^{q)\^ {12) 

Here is the nuclear charge-current operator which connects the initial state | \^o) of the nucleus, 
of energy Eq, with the final states | "^f), of energy Ef, both eigenstates of the {A + l)-body 
Hamiltonian H. The sum runs over the scattering states corresponding to all of the allowed 
asymptotic configurations and includes possible discrete states. G{Ef) is the Green function related 
to H, i.e., 

G(Ef) = . (13) 

Here and in all the equations involving G the limit for 77 — > +0 is understood. It must be performed 
after calculating the matrix elements between normalizable states. 

In Ref. [22] we have showed that the spectral representation of the single particle Green's func- 
tions related to the optical potential allows practical calculations of the hadron tensor components 
for the inclusive (e, e ) scattering. Here, we briefly recall the most important points of the method, 
without discussing in details the approximations involved. Assuming only one-body terms in the 
nuclear current, Eq. (12) can be reduced to a single-particle expression whose self-energy is the 
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Feshbach's optical potential [22]. We then perform a biorthogonal expansion of the full particle- 
hole Green's function in terms of the eigenfunctions the optical potential, V, and its Hermitian 
conjugate, V^ i.e.. 



[£-T-VHE)]\x^i\E)) = () , [E-T-V{E)]\x^,-\e))=0. 
Note that £ and E are not necessarily the same. The spectral representation is 



(14) 



(15) 



and the hadron tensor can be written in an expanded version in terms of the single-partcile wave 
function, | ipn), of the initial state, corresponding to the energy and whose spectral strength is 
Xn, as 



-— Im 

TT ^ 



d£- 



S + ir) 



Ti:^{£,Ef-en) 



(16) 



where 

Ti:^{S,E) = U^r. I f\q)^l-V{E) I x[-\e)){x^s\e) \ ^1 - V{E)f{q) \ . (17) 

In Rcf. [22] we have shown that the factor -s/l — V'{E) accounts for interference effects between 
different channels and allows the replacement of the mean field V with the phenomenological optical 
potential Vl- In Eq. (17) V'{E) is the energy derivative of the mean field potential. After 
calculating the limit for rj ^ +0 Eq. (17) reads 

n 

where V denotes the principal value of the integral. Eq. (18) separately involves the real and 
imaginary parts of Tn^. 

Let us examine the expression of Tn^{£, E) at £ = E = E^ — e„ for a fixed n. This is the most 

important case since it appears in the first term in the right hand side of Eq. (18), which gives 
the main contribution. Disregarding the square root correction, due to interference effects, one 
observes that in Eq. (17) the second matrix element (with the inclusion of -v/A^) is the transition 
amplitude for the single nucleon knockout from a nucleus in the state | ^'o) leaving the residual 
nucleus in the state | n). The attenuation of its strength, mathematically due to the imaginary part 
of V^, is related to the flux lost towards the channels different from n. In the inclusive response 
this attenuation must be compensated by a corresponding gain due to the flux lost, towards the 
channel n, by the other final states asymptotically originated by the channels different from n. In 
the description provided by the spectral representation of Eq. (18), the compensation is performed 
by the first matrix element in the right hand side of Eq. (17), where the imaginary part of V has the 
effect of increasing the strength. Similar considerations can be made, on the purely mathematical 
ground, for the integral of Eq. (18), where the amplitudes involved in Tn^ have no evident physical 
meaning as £ ^ E^— En- We want to stress here that in the Green function approach it is just the 
imaginary part of V which accounts for the redistribution of the strength among different channels. 

The cross sections and the response functions of the inclusive quasielastic electron scattering 
are calculated from the single particle expression of the hadron tensor in Eq. (18). After the 
replacement of the mean field V{E) by the empirical optical model potential Vl(-E), the matrix 
elements of the nuclear current operator in Eq. (17), which represent the main ingredients of the 
calculation, are of the same kind as those giving the RDWIA transition amplitudes of the (e, e'p) 
reaction. Thus, the same treatment can be Tised to describe the initial and final state wave functions 
and the one- nucleon electromagnetic current. 



Rerr (£^/ - £n, Ef - sn) - -V r d£- \ ^Imr„^'^(f , E^ - £„) 

JM Ef - En— O 



,(18) 
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FIG. 4: Longitudinal (upper panel) and transverse (lower panel) response functions for the ^^C(e, e') reaction 
a.t q — 400, 500, and 600 McV/c, respectively. Solid and dotted lines represent the relativistic results with 
and without the inclusion of the factor in Eq. (19), respectively. Dashed lines give the result without the 
integral in Eq. (18). Dot-dashed lines are the contribution of integrated single nucleon knockout only. The 
data are from Ref. [40]. 



A. Results and discussion 

The longitudinal and transverse response functions for ^^C at q = 400 McV/c arc displayed in 
Fig. 4 (left column) and compared with the Saclay data [40]. The low energy transfer values are 
not given because the relativistic optical potentials are not available at low energies. 

The agreement with the data is generally satisfactory for the longitudinal response. The trans- 
verse response is underestimated. This is a systematic result of the calculations and was also found 
in the nonrelativistic approach of Ref. [41]. It may be attributed to physical effects which are not 
considered in the single-particle Green function approach, e.g., meson exchange currents. 

The effect of the integral in Eq. (18) is also displayed. Its contribution is important and essential 
to reproduce the experimental longitudinal response. 

The contribution arising from interference between different channels gives rise to the factor 

^1 - V'l{E) = Vl - PS'{E) - V'{E) . (19) 

We see, however, that here it gives only a slight contribution, due to a compensation between the 
energy derivatives S'{E) and V'{E). 

The contribution from all the integrated single nucleon knockout channels is also drawn in Fig. 
4. It is significantly smaller than the complete calculation. The reduction, which is larger at lower 
values of to, gives an indication of the relevance of inelastic channels. 

The longitudinal and transverse response functions for ^^C at q = 500 (middle column) and 
q = 600 MeV/c (right column) are also displayed and compared with the Saclay data [40] in Fig. 
4. As already found at g = 400 MeV/c, a good agreement with the data is obtained in both cases 
for the longitudinal response, while the transverse response is underestimated. Only a slight effect 
is given by the factor in Eq. (19) arising from the interference between different channels. The role 
of the integral in Eq. (18) decreases increasing the momentum transfer. 

In Fig. 5 wc consider the ^^0(e, e') inclusive cross section data taken at ADONE-Frascati [24] 
with beam energy ranging from 700 to 1200 MeV and a scattering angle — 32°. The NLSH wave 
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FIG. 5: The cross section for the inclusive 1*^0(6,6') reaction at = 32° and Ee = 700, 1080, and 1200 
MeV. The data are from ADONE-Prascati [24]. Line convention as in Fig. 4. 



functions have been used in the calculations. The agreement with data is good in all the situations 
considered. The integral in Eq. (18) produces a reduction which is now essential to reproduce the 
data at 700 MeV, which correspond to a momentum transfer q < 400 MeV/c. Its contribution can 
be neglected when q 2± 600 MeV/c. The effect of the factor in Eq. (19) is very small. 
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